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Introduction 

Seismological calibration of solar models to estimate the age of the Sun necessarily depends 
predominantly on the frequencies of the lowest-degree modes which penetrate into the energy- 
generating core where the greatest evolutionary change in the stratification occurs. Most 
commonly this is accomplished by fitting the asymptotic formula 
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to observed high-order frequencies v an / of order n and degree / to determine the coefficients 
fn, e, and A^j; here L = / + i. The most /—sensitive terms, at each degree 2k — 1, namely 
Ak if, are, on the whole, the most sensitive to core conditions, and the least sensitive to the 
structure of the envelope (cf. Houdek &i Gough 2007b). Therefore it is one or more of these 
that are the best determinants of stellar age. Eq. (1) is valid only if / <C n and n is large, such 
that the spatial scale of variation of the equilibrium state is everywhere much greater than 
the inverse vertical wavenumber of the mode. But that condition is not actually satisfied 
in the Sun: there is small-scale variation associated with ionization of abundant elements 
and the near discontinuity in low derivatives of the density at the base of the convection 
zone, which we call acoustic glitches, and which add the components u gn j to u„j that are in 
general oscillatory with respect to n. Ignoring these components introduces systematic errors 
into a straightforward fitting of Eq. (1) to v an /, errors that are evident in the undulatory age 
estimates as the limits of the frequency range adopted for the fitting are varied (Gough 2001). 
In an attempt to obviate these errors, Houdek &i Gough (2007a, 2008) estimated the glitch 
components u gn j by fitting to second differences (with respect to n) of the observed frequen- 
cies an asymptotic formula designed to represent the base of the convection zone and the two 
ionization zones of helium. In reality there is also an upper-glitch component, produced by the 
ionization of hydrogen and the upper superadiabatic boundary layer of the convection zone, 
which appears to be difficult to model in a reliable manner; when fitting the second differences 
Houdek & Gough (2007a) represented that component, coupled with the second differences 
of Eq. (1), somewhat arbitrarily as a series P(f n /) of inverse powers of v„\. Because the 
upper glitch component is relatively smooth, they subsequently tacitly regarded it as being 
included in the smooth asymptotic expression (1) by adjusting the observed frequencies by 
only the component v gn \. Because the upper glitch is quite close to the surface (partly in the 
evanescent zones of the modes), its influence on the eigenfrequencies is essentially indepen- 
dent of /, and so should not have materially affected the fitted coefficients A k j with j > 0. 

Modification to the calibration procedure 

In the work we report here we have tested the stability of the procedure by including in u s „j a 
representation of the upper-glitch component. To this end we summed the second-difference 
representation P to obtain an estimate of its contribution to the frequencies. There is some 
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Figure 1: Contribution of the upper-glitch 
component to the glitch frequencies v g „j, 
obtained from the BiSON observations u on j 
(Basu etal. 2007), as a function of v on j. 
It was obtained from summation of the 
series P(h> n j), whose coefficients were de- 
termined from fitting to the second differ- 
ences of the observed frequencies the sec- 
ond differences of an asymptotic formula 
representing the glitch components v %n / 
(Houdek & Gough 2008). The upper-glitch 
component is produced by the ionization of 
hydrogen and the superadiabaticity of the 
surface boundary layer. 



ambiguity in how one separates smooth and glitch components near the surface, which is 
exhibited by the two undetermined constants of summation of the second differences; here 
we chose those constants by minimizing the error-weighted sum of the squares of the upper- 
glitch frequencies. The outcome is plotted in Fig. 1 using BiSON data (e.g. Basu etal. 2007) 
up to degree / = 3. 

After fitting Eq. (1), with K = 3, to the resulting glitch-adjusted observed frequencies, 
the coefficients i/q, e and Ak n were found to be naturally somewhat different from the results 
obtained without the upper-glitch adjustment. But the coefficients A^ j, are similar. There is, 
however, a slight difference, which is evidently a product of an inadequacy of the asymptotic 
formulae to reproduce precisely the observed frequencies of the Sun. 

Result 

The result of the present model calibration against BiSON data (e.g. Basu etal. 2007) is 

t Q =4.63 ±0.02 Gy, (2) 
a value in fair accord with our previous estimates (Houdek &c Gough 2007b, 2008). The errors 
quoted here come solely from the stated observed frequency errors, which we have assumed 
to be statistically independent, and take no account of (systematic) errors in our procedure; 
that the value (2) differs from our previous estimates by as much as 2.5 a suggests that such 
systematic errors could be present at a level at least as great as the random errors. Our 
current value for the solar age is lower than the previous estimates by essentially this method, 
although, in contrast to many earlier estimates, it remains greater than the age of Model S 
of Christensen-Dalsgaard et al. (1996), which we used as our reference. It is also greater 
than that of many, if not all, meteorites. We have not yet completed our investigation of the 
robustness of the result, so we offer it still as a preliminary estimate. 
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